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Recent observations by the ROSAT X-ray satellite of the Vela supernova 
remnant 1 have revealed, in addition to the previously identified compact 



nebula 2 , a nearly circular emitting region with a radius of about 4 degrees. 
The Vela pulsar is slightly off the center of this circular region, consistent 
with its measured proper motion 3 of about lOOkrns -1 and an age of about 

. 10 4 yr. The emitting region is bounded by the main supernova shock. Just 

O. 

outside the shock, the X-ray image reveals several well-defined V-shaped 
^ ■ features extending radially outwards 1 . These features are most likely wakes 

Oh; 

produced by objects moving supersonically through the outside medium. 
The shapes and orientations of the wakes suggest that these objects have 
been ejected from the center of the supernova explosion. Their present 

X' 

positions indicate that they have been moving with a mean velocity of a 

a ; 

few thousand kms" 1 . We show that pre-existing objects, such as planets in 
orbit around the progenitor star, could not have been accelerated to suffi- 
ciently high velocities or would have been destroyed. Instead, we propose 
that the observed objects are fragments ejected during the formation of 
the neutron star. Fragmentation during gravitational collapse is a natural 
consequence of both convective 4 and rotational 5 instabilities. 

The X-ray emitting wakes observed by ROSAT appear in the North, East and 
West of the Vela supernova remnant and are stretched radially outwards 1 . Their 



angular distance from the center is about 4°, giving a mean velocity Vf = 2.8 x 
10 3 d,4 r^ 1 kms -1 , where 400 ^pc is the distance to the remnant and 10 4 T4yr is its 
age 3 . The initial ejection velocity must clearly have been > Vf. For the brighter of 
the three features, the ROSAT PSPC countrate is about 0.25 s -1 and the emitting 
area is about 32arcmin 2 , giving an X-ray luminosity L x ~ lO 34 ^ erg s _1 (J. Triimper, 
personal communication). This implies a total radiated energy of ~ 3 x lO 45 ^ 2 ^ erg 
over the lifetime of the remnant. If a moving fragment were to supply this energy 
while decelerating from the above Vf, its mass should be J> 10~ 4 M Q . 

We can immediately rule out the possibility that the fragments are pre-existing 
objects (e.g. planets) that were accelerated by the explosion. The progenitor star 
was most likely a supergiant of radius R p > 10 13 cm. The total energy imparted 
to the ejected stellar envelope in a supernova explosion is E s ~ 10 51 erg. This en- 
ergy includes both thermal and bulk velocity components. The initial velocity of 
the shocked fluid is v s ~ 10 4 kms _1 . Thus, the total momentum flowing out is 
< 2E s /v s ~ 2 x 10 42 gcms _1 . The fraction of this momentum outflow that can be 
intercepted by an object of radius Rf outside the progenitor is < (Rf/R p ) 2 . Objects 
orbiting inside the progenitor are dragged to the center and destroyed on a very short 
timescale 6 ~ 10 3 s. The ejection velocity for an external object with mass Mf and 
mean density pf is Vf < [2E s /(v s Mf)](Rf /R p ) 2 , giving 

c/ < 6km .-. ( - 1 ' 3 (uu-.) - 2/3 f^_) - 2 . (1) 

1 V 10 M o/ \lgcm- 3 / V 10 3 cm/ v; 

For planets this upper limit is at least three orders of magnitude smaller than the value 
ofvf deduced from the observations. Even if the planets could somehow survive inside 
the supergiant envelope, the explosion would destroy them. Indeed, disruption occurs 
whenever the nonuniform pressure forces exerted by the ejected gas exceed a significant 
fraction ~ 0.1 of the internal gravitational force that keeps the planet bound. A planet 
at a distance r away from the center must be accelerated to its ejection velocity > Vf 
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in a time r/v 8 , and the condition Vf/(r/v s ) < O.lGMf /R"j translates to the lower 
bound 

' > 10l5cm {T^kf 3 (l^)" 3 fcnfcr)" < 2 > 

Expressions (1) and (2) show clearly that planets cannot be accelerated to the required 
velocity without being totally disrupted. 

Instead we propose that the observed wakes are associated with fragments that 
were ejected from inside the progenitor core during an asymmetric gravitational col- 
lapse. During the acceleration phase, the density inside the fragments was pi ~ 
10 10-15 gem -3 , high enough for the fragments to avoid disruption. Direct evidence for 
asymmetric core collapse is provided by the large observed proper motions of young 
pulsars . Two specific mechanisms can result in the formation of small fragments. 
Rotationally-induced instabilities during gravitational collapse 5 can lead to the devel- 
opment of either an axisymmetric, self-gravitating disk 8 , or an ellipsoidal deformation, 
which then leads to mass shedding through outgoing spiral arms 9 ' 10 . The spiral arms 
can later fragment through a sausage instability, leaving small compact objects or- 
biting the central neutron star core 11 . The fragments in this case would be initially 
at close to nuclear density, pi ~ 10 14 gcm -3 , and ejected from ri ~ 10 1_2 km. They 
would be formed just outside the radius where the bounce shock forms initially, and so 
they could be shock-accelerated and ejected before the shock stalls 12 . Alternatively, 
Rayleigh- Taylor (convective) instabilities in the outer iron core 4 can lead to small blobs 
of overdense material forming behind the shock front. These overdense blobs could 
be accelerated by absorbing a small fraction of the neutrino flux, i.e., through the 
same mechanism that is now commonly thought to re-energize the stalled shock and 
power the ejection of the stellar envelope 4 ' 12 . In this case the initial density is lower, 
Pi ~ 10 10 gcm -3 , and the ejection is from ~ 10 2-3 km. 
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In response to the release of nuclear energy in their interior, fragments that are 
too small may disintegrate. Consider a fragment made of hot (T J> 1 MeV) , highly neu- 
tronized material. As the fragment expands, nucleons recombine quickly into nuclei. 
At relatively late times the energy release is dominated by /3-decays 13 . The maxi- 
mum amount of nuclear energy released is about 8 MeV per nucleon if all the matter 
transforms into iron. The actual energy which gets thermalized inside the fragment is 
likely to be considerably smaller, because part of the nuclear energy is carried away by 
neutrinos, and because the fragment may not form out of highly neutronized matter. 
Ignoring degeneracy pressure, a necessary condition for a fragment to remain gravita- 
tionally bound is that the thermalized nuclear energy be smaller than the gravitational 
binding energy per baryon, 

We therefore obtain a conservative estimate of ~ 0.01M© for the minimum mass of a 
fragment that would not disintegrate as a result of nuclear energy release. 

Even if it can sustain the nuclear energy release, a fragment could still be unbound 
if it is produced with an initial temperature much larger than the virial temperature. 
Just like a nascent neutron star, a newly formed hot fragment can cool rapidly by 
neutrino emission, to which it is optically thin. When the neutrino cooling time 
becomes comparable to the hydrodynamic time, the fragment begins expanding and 
may become unbound if its temperature is still too high. The neutrino cooling time 
due to pair annihilation is 14 10~ 4 s (T/10 MeV) 5 . If we take the sound speed to be 
~ 0.1c and the initial fragment size to be ~ 1km, then the hydrodynamic expansion 
time is ~ 3 x 10~ 5 s. This becomes comparable to the neutrino cooling time at T = 

1 /5 

T c ~ 12 MeV. The exact value of the sound speed is not important since T c oc c s . 
Using equation (8) we conclude that for initially hot fragments, the mass must be 
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O.O1M to avoid disintegration. This estimate turns out to be comparable to that 
obtained above by considering only the nuclear energy release. 

As the fragments move away from the center of the explosion on ballistic orbits, 
they expand and eventually settle into a gravitational equilibrium with a density ~ 
lgcm -3 . However, they could remain hot for a time ^> 10 4 yr because of their large 
initial thermal heat content as well as the energy release from radioactive nuclei in 
their interiors. An evaporation process should result. We write the evaporation rate 
in terms of the mass flux at the surface, 

-M f = 4np e R 2 f v e , (10) 

where p e is the density of the wind at ejection and where the ejection speed v e must 
be larger than the escape velocity from the surface of the fragment, 

«. Z >W S 1.3 x 10* km " 3 ( Ti ^ FJ ) . (11) 



From mass conservation the density in the terminal velocity outflow is 



peR 2 f 

P(r) « (12) 



We can estimate the radius where the outflow will be deflected by the ram pressure 
of the external medium from the condition 

P{vt ~ Vesc) = PextVf, (13) 

where p ext ~ 10~ 24 g cm -3 is the ambient interstellar density. Taking (i>g— i>g SC ) ~ v 2 sc , 
we find 

/ \M f \v \ lf \ 4x 16 / _A^_A V2 ( vj_ \ - 1 

\^ Pext v)) \M max ) VSxlOSkms" 1 ; 1 } 

where — M max = lO _6 M0yr _1 is the maximum evaporation rate allowed for a frag- 
ment of mass 1O _2 M over the current pulsar lifetime. This radius ra gives the effective 
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cross-section for estimating the drag force F d ~ Pext^j x (7rr^). The corresponding 
energy dissipation rate F^vj yields a luminosity, 

L x « 7rp ex ^?r3 ~ 4 x lO^ergs" 1 ( ^ ] [ \ . (15) 

fd \2 x 10 3 kms -1 / yM max y 1 J 

This crude estimate comes very close to the observed X-ray luminosity determined from 
the ROSAT observation. Most of the energy is expected to be released by thermal 
bremsstrahlung in X-rays since the post-shock temperature of the gas is high. In fact, 
the Vela supernova may have ejected additional high velocity fragments. However, 
fragments with \M/\ <C \M m ax\ are too faint to be detectable by ROSAT, while those 
with a much higher evaporation rate would have disappeared by now. 

To power an evaporation rate close to M max requires little energy deposition at 
the surface of a fragment, namely ~ 100 eV per baryon. This energy can be easily 
supplied over 10 4 years by short-lived radioactivity or by the initial heat content of 
the fragment. In order to maintain an appreciable evaporation rate, the temperature 
at the surface of the fragment must be ~ 1/3 of the escape temperature 15 , 

T s ,3x 10 5 K ( -^-) 2/3 ( • (16) 

Vio - 2 m q J Vigcm- 3 ; v ; 

The material in the evaporative outflow is initially fully ionized. As the surface material 
expands and rarefies to become part of the supersonic wind, it cools adiabatically with 
T oc p 2 / 3 . The gas starts to recombine within a distance ~ Rf and eventually reaches 
a sufficiently low ionization level, so that the outflow becomes optically-thin, at a 
temperature T p h ~ 5000K around r p h ~ 3i?/. The emerging optical luminosity from 
the wind photosphere is then, 

where a is the Stefan-Boltzmann constant. This estimate is highly uncertain because 
of the unknown contribution of heavy elements and molecular absorption to the opacity 
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of the wind. Nevertheless, it shows that the optical luminosity of the fragments may 
be detectable. It should therefore be useful to search for optical emission near the 
leading edges of the X-ray emitting wakes. 
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